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Abstract. The center of the Milky Way hosts a massive black hole. The 
observational evidence for its existence is overwhelming. The compact radio 
source Sgr A* has been associated with a black hole since its discovery. In the last 
decade, high-resolution, near-infrared measurements of individual stellar orbits 
in the innermost region of the Galactic Center have shown that at the position 
of Sgr A* a highly concentrated mass of 4 x 10 6 M is located. Assuming that 
general relativity is correct, the conclusion that Sgr A* is a massive black hole 
is inevitable. Without doubt this is the most important application of stellar 
orbits in the Galactic Center. Here, we discuss the possibilities going beyond the 
mass measurement offered by monitoring these orbits. They are an extremely 
useful tool for many scientific questions, such as a geometric distance estimate to 
the Galactic Center or the puzzle, how these stars reached their current orbits. 
Future improvements in the instrumentation will open up the route to testing 
relativist ic effects in the gravitational potential of the black hole, allowing to 
take full advantage of this unique laboratory for celestial mechanics. 



1. Observations 



Observationally, the Galactic Center (GC) region is characterized by the high 
extinction screen along the line of sight through the Milky Way plane. In order to 
observe stars in the GC, one needs to go to the near-infrared (NIR). Given the 
small angular scales and the crowding of stars, high resolution methods have 
been the clue for detecting the stellar orbits. In the 1990s NIR observations 



were using the Speckle technique (Eckart & Genzel 1996; Ghez et al. 1998). 



Since 2002 adaptive optics observations, both for imaging and spectros copy are 
performed r outinely, either using natural guide stars or laser guide stars (Schodel 
et al.||2002| |Eisenhauer et aL||2005||Ghez et aL||2005afl . The images obtained at 
the diffraction limit of the telescope (#fwhm ~ 55 mas for the K-band around 
2.2 /im at an 8 m-telescope) reveal w 100 stars within r = 1" of the massive 
black hole (MBH, figure [I]). The brightest of these so called S-stars, S2, has 
an apparent magnitude of mx = 14. The high density of stars limits how 
faint a star one can detect in the central arcsecond. Under good conditions, 
stars down to mx = 18.5 can be detected. The astr ometric precision from 
such data is in the 200 — 300 ^as range for bright stars (Ghez et al. 2008; Fritz 



et al. 2010). Spectroscopically, the use of adaptive-optics assisted integral field 



1 



2 Gillessen et al. 




Figure 1. Finding chart of the innermost region of the nuclear star cluster 
in the Galactic Center. The image is based on a natural guide star adaptive 
optics image from NACO at UT4 (Yepun) of the VLT on 2007 July 20 in 
the H band. Stars that are unambiguously identified in several images have 
designated names, ranging from SI to SI 12. Blue labels indicate early- type 
stars, red labels late-type stars. Stars with unknown spectral type are labelled 
in black. 



spectroscopy is state of the art. In the central arcsecond, late-type stars as 
faint as ttik ~ 17 can be identfled, given the prominent CO absorption features 
around 2.3 /im. For early-type stars, the current limit is tuk ~ 16.5, and the 
detection of the Brackett-7 line is hampered by the presence of line emission 
from the background gas. Radial velocities can be measured with a precision 
of w 15km/s for the brightest early-type stars and for many late-type stars, 
limited by the accuracy of the wavelength calibration. For fainter early-type 
stars, the errors are dominated by the signal-to-noise ratio and typical errors are 
between 50 — lOOkm/s. From these a strometric and spectroscopic measurements 
the orbits of ~ 30 star s are known flGhez et al.||2005b ; Eisenhauer et al. 2005; 
Gillessen et al.|[2009b| ). Surprisingly, a large fraction of these stars deep in 



the potential of the MBH ar e young. Their ex istence in the place where they 
are observed poses a puzzle (Ghez et al. 2003). Due to the tidal forces of the 
MBH they cannot have formed in-situ, but their age is much smaller than the 
migration time scale. 
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Figure 2. The orbit of S2. Left: |Ghez et al.| (|2008) (red, open symbols) and 



Gillessen et al. ( 2QQ9b[ ) (blue, filled symbols) have published orbital data for 
S2. Simply overplotting the two data sets and orbit fits reveals an apparent 
mismatch larger than t he precision of the individual measurements. Middle: 
jGillessen et al.| ( |2QQ9a| showed that the difference can be understood com- 
pletely by the uncertainties of the two coordinate systems. A simple offset 
in position and velocity is sufficient to make the two data sets agree. Right: 
The radial velocity data. 



The orbit of S2 



Since the beginning of the observations in 1992, S2 has completed a full revo- 
lution around Sgr A* (figure [2]). Its orbital period is 15.9 years, the sem i-major 
axis of the orbit is 125 mas dGhez et al.|[2008l |Gillessen et~aT1|2009b|a| ). omce 
the proper motion is measured in angular units per time, but the radial velocity 
in physical length units per time, an orbit model contains the distance to the 
system as parameter. Hence, S2 allows for a geometric determination of i?o, 
the distance to the GC (Salim & Gould 1999). The current best such estimate 
(IGillessen et al.||2009a|) is 



Ro = 8.28 ± 0.15| sta t ± 0.30| sys kpc . 



(1) 



Its accuracy is limited by systematic uncertainties. The corresponding value for 
the central mass is 



M = 4.30 ± 0.06| stat ± 0.35| Ro x 10 6 M 



(2) 



The first error gives the uncertainty if the value of Ro were known perfectly 
well, the second, dominant term is the error induced by the uncertainties on Rq. 
This indicates that the M and Ro are strongly correlated, empirically we find 
M oc Rq. The exponent here would be 3 for a purely astrometric data set and 
for a purely spectroscopic one. 

The systematic unc ertainties had been underestimated in previous works 
(Eisenhauer et al. 2003). A large part of them comes from the fact that it is 
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hard to es tablish an absolute reference frame in the NIR in the GC region (Reid 
et al.||2007 ). This is illustrated by the fact, that the published data sets do not 
agree at first glance (figure [2]). However, allowing for an offset in position and 
velocity betwee n the two coordinate sy stems suffices to bring the measurements 
into agreement ( |Gillessen et al.||2009a ). 

The NIR coordinate system is relying on the astrometric measurements of 



SiO maser stars that are visible both in the NIR and in the radio regime (Reid 
et al. 2007). By construction, the radio source should be at rest at the origin 
of the NIR coordinate system. The position of the central mass observationally 
is a distinct quantity and thus it is worth to compare how well the positions of 
mas s and Sgr A* do agree. We find that they c oinc ide to within less than 3 mas. 
~ (12005b, 
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Ghez et al.|fl2008D and |Gillessen et aL|fl2009b| ) have 
e the S2 orbit can be in an extended configuration. 

5% of the mass of the MBH can 

For 



tested, how much mass insid< 
These studies consistently find that at most 

reside at radii between the peri- and the apocenter distance of the S2 orbit, 
example, Gillessen et al. (2009b) find that the fraction 77 of the extended mass 
to the MBH mass is 



V = 0.018 ± 0.014| sta t ± 0.005| mode i 



(3) 



The measurement uncertainty is still too large to give astrophysically rele- 
vant constraints. Various arguments point towards a value of 77 ~ 10 -3 — 10 -4 . 

• The drain limit gives the maximum number of objects that in a steady 
state situation can reside close to a MBH. From the work of Alexander 



& Livio (2004J) follows that if the dark mass consists of stellar black holes 
with 10 M©, one would expect 77 < 1.1 x 10 -3 



• Theoretical modelling, either in a Fokker-Planck approximati on or using 
explicit N-body simulations yield c onsistently 77 « 5 x 10 -4 (Hopman & 
Alexander||2006l |Freitag et aL||2006| ) . 



The overdensity of X-ray binaries in the central parsec ( Muno et al.||2005 ) 
also estimates the number of stellar black holes in the GC. The m easured 
density corresponds to 77 « 2 x 10 -4 ( Deegan fc Nayakshin||2007 ). 



The amount of cosmological dark matter is completely negligible compared to 
the expected astrophysical backgroun d density. The dark matte r density is at 
least two order of magnitudes smaller ( Vasiliev fc Zelnikov||2008 ) . 

The only hypothesis which can be tested in the near future is, whether the 
S-stars have reached their orbits via normal two-body relaxation from a birth 



place in the region of the disks of young, massive stars (Paumard et al. 2006 



Lu et al.||2009[ jBartko et al.||2009| ). This would require significantly more mass 
than what is seen in stars. Stellar black holes could lower the relaxation time 
such that the lifetime of a B-star would suffice to migrate from the disk region 
into the central arcsecond. The corresponding value of 77 is « 0.03, relatively 
close to what can be measured. However, such a scenario is disfavored anyhow. 
The prospects for detecting relativistic effects are much better. |Zucker et al 
showed that the relativistic effects of order /3 2 in the radial velocity are 



fl2006| ) 



detectable in the orbit of S2 with today's technique during a pericenter passage. 
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Figure 3. Contributions of various effects to the errors of astrometric po- 
sitions as a function of brightness of the star considered. 



Besides the classical Romer effect, these are the transverse Doppler effect and 
gravitational redshift. Here, it is important to realize that it is not sufficient to 
compare for a given fixed orbit the Keplerian and relativistic fits. Instead, one 
needs to compare the two fits for a given set of data, since the orbital elements 
need to be determined from the same data. The relativistic effects are degenerate 
to a strong degree with the orbital elements when fitting orbits. 

The largest effect detectable in astrometry is the prograde pericenter pre- 
cession due to the Schwarzschild correction of the Newtonian potential. For 
S2, the difference between best fitting Keplerian and relativistic orbit is smaller 
than 0.5 mas at all times for the next years, and it will be considerably more 
difficult to detect the precession than the effects on the radial velocity. This can 
also be seen from precision by which the position angle of the ascending node 
is determined. Currently, after 18 years of monitoring, Acj — 0.84°, while the 
precession amounts to 0.22° per revolution (|Rubilar fc Eckart||2001[). 



3. The limits of astrometry 



Fritz et al. (2010) have investigated the limits of adaptive-optics based astrom- 
etry in the crowded GC region. The results are summarized in figure [3j For 
the brightest S-stars, the limitation currently is the knowledge of the optical 
distortions of the camera. For fainter stars, the uncertainty is dominated by the 
so-called halo noise. The term refers to the fact that the imperfectly determined 
seeing halos of the brighter stars perturb the positions of the fainter stars. The 
effect depends on the level of crowding (and thus on distance to Sgr A*) and 
most importantly on the Strehl ratio. The signal-to-noise ratio is not the limiting 
factor. Obviously, higher Strehl ratios would improve the astrometric precision, 
and the same is true if one would be able to increase the spatial resolution. 
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Figure 4. Illustration of the stellar orbits determined so far. Here we show 
the innermost 20 of the known orbits. 



The next big step in angular resolution is expected to come from NIR 
interferometry. The resolution which the VLT interferometer can achieve is as 
good as 3 mas in the the K-band. Due to the baselines of « 100 m length, 
this exceeds even the performance of the next generation of extremely large 
telescopes. At that resolution, one should be able to observe stars within the 
central 100 mas that would have orbital periods of w lyear. The relativistic 
precession would amount to several degrees per revolution and should be easily 
detectable within few years of monito ring. The number of stars expected from 
extrapolating the luminosity function ( Bartko et al.|20T0 ) and the radial density 
profile (IGenzel et al. 1120031 ISchodel et al.||2007p in that radial range is between 
three and ten for magnitudes between mx = 17 and 19. 



4. Statistical properties of orbital elements 

Given that w 30 orbits have been determined so far (figure B , one c an start to 
look into the distributions of orbital elements (Gillessen et al. 2009b). They are 
extremely valuable clues for the enigmatic formation of the S-stars. 

The most striking feature in the distribution of orbital elements is the dis- 
tribution of orbital planes. While for the central S-stars the distribution is 
consistent with a random distribution, we find a group of six orbits at larger 
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Figure 5. Orientations of the orbital planes for the known orbits in the GC. 
The blue dots represent stars that belong to the clockwise rotating stellar disk. 
The red dots are the other stars, which appear to be randomly distributed. 



radii (r « 1") that share a co mmon plane. This explicitly confirms the existence 
of the clockwise stellar disk ( Paumard et al.||2006 ; Lu et al. 2009; Bartko et al 
|2009[ ) that previously was defined only in a statistical sense. The six orbits all 
have a moderate eccentricity (e w 0.4) and the brightness of these stars is larger 
than for the stars at smaller radii. We indeed see two populations, the disk 
stars (more massive, orbiting in a common plane, moderate eccentricity) and 
the B-type S-stars (randomly oriented orbits with higher eccentricity). 

Secondly, the distribution of eccentricities (of non-disk stars) is slightly 
hotter than for a thermal (n(e) oc e) distribution: 



n(e) oc e 



2.8±0.8 



(4) 



This potentially is a key to the formation of the S-stars. In the Hills mechanism 
pillsl[l99l| ) scenario, massive perturbers (for example molec ular clouds) scat- 
ter field stars into a near loss-cone orbit (Perets et al. 2007). If that happens 



to a binary, the three-body interaction with the MBH at closest approach can 
break the binary, leaving one member as S-star and ejecting the other star at 
high velocity, which might explain the existence of hypervelocity stars ( |Brown| 



et al. 2005, 2006). Such a scenario indeed would yield an ec centricity distribu 



tion hotter than thermal, however it would be much hotter (Ginsburg &; Loeb 



2006, Alexander 2007). Hence an additional relaxation mechanism needs to 



be invoked, maybe resonant relaxation (Hopman &; Alexander 2006). Clearly, 
determining the eccentricity distribution more accurately is important to the 
binary-capture scenario. 

The distribut ion of semi-ma j or a xes tests the density profile: n(r) oc r a 
for n(a) oc a a+2 (Schodel et al. 2003). Our measurements for r < 0.5" yield 
n(r) oc r -1,1 , which agree s well with previous works at larger radii (Genzel et al. 



20031 ISchodel et al.|[2007| ). 
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Figure 6. . Energy dependence of the absorption of the TeV radiation from 
Sgr A* due to a star with given tem perature for an angular separation of 
1 mas. From Abramowski et al. (2009). 



5. Constraining the presence of an intermediate mass black hole 



An intermediate mass black hol e (IMBH) in the GC is proposed as alternative t o 
the binary capture mechanism ( jMerritt et al.||2009 ; Gualandris &; Merritt 2009). 
Also, an IMBH might play an important role for the young, massive stars. If 
they have reached the GC as part of an inspiraling cluster, the IMBH might 
have helped to bind the cluster such that it was able reach dista nces of a few 
arcseconds. (Hansen & Milosavljevic 2003; Giirkan fc Rasio||2005 ) . 

The current limits on the motion of Sgr A* obtained from stellar orbits are 



not competitive with the radio measurements of the motion of Sgr A* (Reid &; 



|Brunthaler||2004| |Merritt et all [2009) since they only exclude a region in the 



parameter space of mass and distance of the IMBH at relatively small distances 
and large masses that is excluded form the radio data already. 

( |2008[ ) showed than an IMBH would change the orbital ele- 
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ments of the S- stars continuously, and maybe by an amount which is not much 
below the current measurement accuracy. Hence, testing the orbital elements 
for constancy can probe the presence of an IMBH. In a more extreme form, one 
might also hope to be lucky to observe a gravitational encounter between an 
IMBH and an S-star. 



6. The nature of the TeV source in the Galactic Center 



At th e position of Sgr A* also a TeV source is observed ( Aharonian fc et al. 

2004) . It is unlikely that the emission is due to dark matter annihilation (Horns 

2005) . Instead, it probably is associated either with Sgr A* itself or a pulsar 
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wind nebula a few arcseconds north-west of Sgr A* (Acero &; et al. 2009). The 
pointing uncertainties of the Cerenkov telecopes make it unlikely that from the 
direct measurements one will be able to decide which source is the origin for the 
TeV signal from the GC. 

The S-stars offer a neat probability to test, whether Sgr A* is the TeV 
source. The (UV) radiation field around a star is opaque to TeV gamma rays. 
Hence, if an S-star cross es the line of sight to Sg r A*, one might observe a 
decrease of the TeV flux. Abramowski et al. ( |2009| ) have calculated the size of 
the effect for the known S-stars (figure [6]). A flux decrease of a few percent over 
the course of a few months would be the typical signal of such an pair production 
eclipse. Future TeV detectors should be able to detect this. The absence of the 
signal in turn would exclude that Sgr A* is the TeV source. 
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